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ABSTRACT 

The central compact objects are a newly-emerging class of young neutron stars near the 
centre of supernova remnants. From X-ray timing and spectral measurements, their 
magnetic fields are determined to be ^ 10^" — 10^^ G, which is significantly lower 
than that found on most pulsars. Using the latest electrical and thermal conductivity 
calculations, we solve the induction equation to determine the evolution of a buried 
crustal or core magnetic field. We apply this model of a buried field to explain the 
youth and low observed magnetic field of the central compact objects. We obtain 
constraints on their birth magnetic field and depth of submergence (or accreted mass). 
Measurement of a change in the observed magnetic field strength would discriminate 
between the crustal and core fields and could yield uniquely the birth magnetic field 
and submergence depth. If we consider the central compact objects as a single neutron 
star viewed at different epochs, then we constrain the magnetic field at birth to be 
~ (6 — 9) X 10^^ G. A buried magnetic field can also explain their location in an 
underpopulated region of the spin period-period derivative plane for pulsars. 
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1 INTRODUCTION 

The magnetic fields on the surface of neutron stars (NSs) 
are known to span a large range: from B ~ 10* — 10® G 
for millisecond pulsars and NSs in low-mass X-ray binaries, 
through lO^'^ — 10^'^ G for normal radio pulsars, to 10^** — 
lO'^^ G for magnetars. The primary method for determining 
the magnetic field is by measuring the spin period P and 
spin period derivative P. Then, by assuming that the pulsar 
rotational energy is lost through magnetic dipole radiation, 
one obtains the relation for the magnetic field at the pole 



pp = IbI 



^ Bp = 6.4 X 10'® G (PP) 



1/2 



where 
10-*" i 



7 



Re I A 



'sin^ a/Sc"/ 
R and I are the 



(1) 

4.884 X 
NS radius 



and moment of inertia, respectively, a is the angle between 
the NS rotation and mag netic axes, and Pe = P /10^ cm 
and J45 = J/10*'^ R cm^ (iGunn fc Ostriked [l969l : see also 
IShapiro fc Teukolskvl Il983l : ISpitkovskvll2006l '). Note that a 
coeflicient of 3.2 is often used in the literature, so that the 
magnetic field that is inferred is the equatorial value; since 



email: wynnho@slac.stanford.edu 



we are modelling the field evolution at the pole, we will 
hereafter only refer to the magnetic field at the pole. 

The evolution of NS magnetic fields is uncertain. Radio 
pulsars are seen to span a large range of ages (~ 10^ — 10® yr), 
magnetars are relatively young (~ lO'' yr), and millisecond 
pulsars are old (~ 10* — 10® yr). The diffusion and decay of 
the magnetic field occurs on the Ohmic timescale 

-rohm = A-KoL^/c ~ 4x10^ yr (a/lO^'' s^') [L/l kmf , (2) 

where a is the electrical conductivity and L is the lengthscale 
over which the decay occurs; here we have taken L to be ap- 
proximately the size of the NS crust (see Fig. [T|. Note how- 
ever that there is a large range in rohm due to the density de- 
pendence of the conductivity and the large density gradient 
in the crust. Thus NSs born with B > 10^^ G survive for long 
times at these field strengths. A manifestation of these high 
fields is the observed pulsar/magnetosphere activity seen in 
young (and old) NSs, such as in the Crab pulsar which is 
only ~ 10"^ yr old. The low magnetic fields which exist on 
millisecond pulsars (and on the NS primaries in low-mass 
X-ray binaries) are thought to be due to burial of the field 
by mass accretion from a companion star. However, several 
young (< 10* yr) NSs are believed to have B ~ 10'°-10" G; 
these NSs are the so-called central compact objects (CCOs), 
i.e., NSs near the centr e of their respective supernova rem- 
nants (SNRs; see, e.g., iDe Lucall2008l : iGotthelf fc HalpernI 
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I2OO8I . for review; see also iHalpern fc Gotthelj I2OIOI ). In 
particular, X-ray timing measurements reveal pulsations 
and spin-down rates that imply from eq. ([1} that Bp < 
4.0 X 10^^ G for P SR J0821-4300 in SNR Puppis A 
iGotthelf et al.ll2010l ). Bp < 6.6 x lO" G for IE 1207.4-5209 
in SNR PKS 1209-51/5 2 (also known as G296.5-f 10.0; 
iGotthelf fc HalpernI [2OO7I ), and Bp = 6.1 x 10^° G for 
PSR J 1852-^0040 in SNR Kes 79 (|Halpern fc Gotthelj 
l2010h : hereafter, we refer to these three CCOs as Puppis A, 
IE 1207, and Kes 79, respectively. In addition to timing 
measurements, the X-ray spectrum of IE 1207 shows ab- 
sorption features which could be due t o electron cyclotron 
resonance at B ^ 7 - 8 x 10^° G (jSanwal et al.l |2002| : 
IMereghetti et al.ll2002bl ). while a possible line in the spec- 
trum of PSR J0821-4 300 could indic ate 5 8- 10 x 10^° G 
jGotthelf fc Halpern| [2b09: Sulcima nov et aLlbOld ). Finally, 



and l2.4|l and parameters that are particular to the properties 
of the CCOs. By comparing the measured properties of spe- 
cific CCOs with our calculations, we obtain constraints on 
the birth magnetic field and submergence depth in each NS. 
The commonality of these two values in the various sources 
further strengthens the distinctiveness of this class of NSs. 

In Section (2] we discuss the magnetic field evolution 
equation and the input physics, in particular, the conductiv- 
ities, equation of state, and temperature evolution. In Sec- 
tion [S] we describe our general results. In Section [l] we ap- 
ply our results to the three CCOs with known spin periods 
and period derivatives or magnetic fields, in order to deter- 
mine the strength of their birth magnetic fields and depth 
of submergence. We summarize our results and discuss their 
implications in Section [S] 



ber of this class, the NS in SNR Cassiopeia A (IMurrav et al. 


2002 


: IMereghetti et al.ll2002al: Ransoml2002l:|Pavlov fc Luna 


200£ 


: Halpern fc Gottheli 201Q). X-rav spectral fits sug- 


gest B < 10" G ( 


Ho fc Heinkd l200g|: iHeinke fc Hdl2010l: 


Shternin et al.ll2011 


). 



In this work, we show that young 

(< 10* yr) NSs, 

like the CCOs, can be understood as NSs which are born 
with B > 10^^ G but whose fields have been buried deep 
beneath the surface jRomanil Il990l l. perh aps by a post- 
supernova episode of hypercritical accret ion (jChevalieill 19891 : 
iGeppert et"al] Il999l : iBernal et al]|2010l '). These fields then 
diffuse to the surface on the timescale of 10^ — 10'' yr, so 
that only now do we see a surface B ~ 10^° - lO" G. It 
is important to emphasize that we are primarily concerned 
with field growth in < 10* yr, which is much shorter than the 
estima te given by eq. ([21). We briefly summarize relevant past 
works. lYoung fc ChanmugamI (|l995h studied field growth in 
old NSs (with low temperatures) with buried magnetic fields 
and found that the surface field is restor ed after ~ 10^° yr 
if the total accreted mass is < O.O4M0. iMusUmov fc Pagel 
(|l995h examined field evolution when the field (confined 
to the crust) is buried at shallow depths, where the den- 
sity is psub ~ 5 X 10^'^ and 4.3 x 10^^ g cm^^; as a result 
of the shallow submerg ence, the diffusion times are short 
(a few X 10^ - lO''^ yr). iMuslimov fc Pagel (|l996l ) modelled 
field evolution in three young (< 1900 yr), fast rotating {P = 
30, 50, and 150 ms and initial spin period Pi = 20 — 40 ms) 
radio pulsars with a measured braking index [see eq. (|22|l ]. 
taking into account cooling of the NS, which influences the 
electrical conductivity and hence the field diffusion timescale 
[see eq. ([5}]. The magnetic flux is assumed to be frozen into 
the core on timescales < 10® yr. Only shallow submergence 
[psub ~ (1 — 3) X 10'° g cm~^] is considered since the pul- 
sars require rapid field growth (to B ~ lO'^^ G) given their 
young age; note that the rapid field growth also results in 
a more rapid spin-down rate. Finally. ICeppert et al.l (|l999l ) 
considered field growth of a purely crustal field at three sub- 
mergence depths, Psub = 10^^, 10^^, 10" g cm"^. 

Discoveries (especially the recognition of the CCO class 
of NSs) since the works discussed above, as well as advance- 
ments in the theory of the relevant electrical and thermal 
conductivities, motivate the current work. We follow a sim- 
ilar methodology as in the previous works tha t examined 
magnetic field growth (see ICeppert et al. I I 19991 . and refer- 
ences therein), but here we use updated physics (see Sec. 12.21 



2 NEUTRON STAR MODEL 
2.1 Magnetic field evolution 

Our calculation of magn etic field evolution is based on 
lUrpin fc Muslimovl (|l992h . To determine the evolution of 
the submerged field, we solve the induction equation 



— = -V X - — V X B 

dt \4TTa 



-V X (r;V X B) , 



(3) 



where rj = c? /Ana is the magnetic diffusivity. Our primary 
interest is in the NS crust, which is predominantly in a solid 
state (see below); therefore we neglect internal fiuid motion. 
We assume the magnetic field in the stellar interior is dipolar 
and given by the simple azimuthal vector potential, A = 
BtR^ s{r,t) sm9/r (j), so that the axially symmetric field is 
given by 



2s 

B,R-" ( — cos^r- - ^ 



19s . , 
— — sm t 

r or 



— B-t (^br COS Of + bo sin ( 



(4) 



where _B, is the "birth" magnetic field and the second equal- 
ity defines the normalized fields 6^ and bg. The birth mag- 
netic field is understood to be the surface magnetic field 
after NS formation but prior to mass accretion; accretion 
then buries and compresses the birth field. Assuming only 
radial variations of the electrical conductivity, i.e., a = a-{r), 
eq. reduces to the one-dimensional equation 



ds 



o s 

Qj.2 



2s 



(5) 



and the field evolution is solely determined by the electrical 
conductivity. 

We assume that the accreted material is non-magnetic, 
so that the initial field decreases very rapidly at shallower 
depths than psub- A pre-existing field at these depths would 
shorten the growth time; however this field must be ordered, 
whereas (hyperc r itical) accretion is likely to be turbulen t 
(jChevalieil Il989l : ICeppert et aD Il999l : iBernal et all I2OI0I I. 
The boundary conditions are then 



ds s 
Or R 



(6) 



at the surface r = R and s — >■ constant in the deep interior. 
The constant determines the two internal configurations for 
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the initial magnetic field that we consider (see thick lines 
in Fig. [2]): (1) confined crustal field, i.e., the diffusion of 
the magnetic field does not extend into the core (e.g., as a 
result of a superconducting core) and the field is negligible 
throughout the star except at the submergence density p^ub 
or position rsub, where 



s(r,0) = (i/2)e-[('-''-'')''*'-l' 



(7) 



and Sr = 10 ^rgub; (2) deep crustal "core" field, i.e., the 
field is finite at p > psub or 



s(r,0) = 



1 



for p < psub 
for p > psub 



(8) 



The Ohmic diffusion timescale [see eq. ([2])] near the core 
is much longer than the times under consideration here 
{t < 10* yr), so that a constant core field is valid; also field 
decay and flux expulsion due to superfluid motion in the 
core only occurs at t > 10* yr s ince the low surface magneti c 
field produces little spin-down (|Konenkov fc Geppert|[200ll ). 
Note that the factor of 1/2 in eqs. (O and (O has been in- 
cluded so that the normalized magnetic field is 6r(Psub) ~ 1 
and B(psub) ~ B, [see eq. ©]. 



2.2 Electrical conductivity 

The dominant contribution to the electrical conductivity de- 
pends on whether the matter is in a liquid or solid state. The 
melting temperature Tm is given by 

= 3.04 X lO'' K (Z/26)^''^(170/r)x, (9) 

where x = {Zpfi/AY^^, Z and A are the charge number 
and mass number of the ion, respectively, T = Z e /akaT, 
a = (3/47rni)'^'''^, ni is the ion number density, and pg = 
p/10^ g cm"'^. The electrical conductivity is mainly deter- 
mined by electron-ion (ei) scattering when T > Tm and 
electron-phonon (ph) scattering when T < Tm. For illus- 
trative purposes, these are given by 



8.53 X 10^^ s"^ 



1.21 X 10'''' s 



^Aei(l-fa;2) 
.1 Vu'^ + 0.017 



Tu 



2 + 



(10) 
(11) 



where A^i is the Coulomb logarithm, u = 0.45T/Td, and Td 
is the Debye temperature 

Td = 3.4 X 10*^ K {Z/A)^/\^^\ (12) 

Equations (fTOl) and (fTT)l (|Yakovlev fc UrpinI [l980l : see also 
lltoh et all I1993I . and references therein) were used by all 
previous works to solve the induction equation in order to 
study magnetic field growth (see Sec. [!}. Instead, we use 
here CONDUCTO^H) which implements the latest advance- 
ments in calculating th e conductivities , including electron- 
electron contributions (*Potckhi n et al.l 1 19991 : ICassisi et al.l 
I2OO7I : IChugunov fc Hacn scl 200 3); we also test ed SPIT- 
TING (N. Itoh, private comm.; lltoh et all 120081 '). We as- 
sume no contribution by impurity scattering, which only 
becomes important at high densities and low temperatures, 
and no magnetic field effects on the conductivities. At 



http: / /www. ioffc.ru/astro/conduct/ 



high magnetic fields and low temperatures, electron mo- 
tion is strongly influenced b y the magnetic field (see, e.g., 
lYakovlev fc KaminkeJl 19941 ). As a result, the electrical and 
thermal conductivities become anisotropic, depending on 
whether the motion is paral lel or transverse to the direction 
of the field (|Potekhinlll999l ). We only consider here young 
NSs with high interior temperatures; thus magnetic field ef- 
fects are minimal since the field is in the non- or weakly- 
quantizing regime, i.e., T > Tb, where 



Tb = 1.34 X 10** K (B/10^^ G)(l + x^) 



-1/2 



(13) 



and p > pB = 2.228 x 10^ g cmr'\A/ Z){B /IQ^''^ Gf^"^. 
However, magnetic field effects can be important in the 
cooler, low-density atmosphere, which pro duces the ob- 
served X-ray emission from NSs (see, e.g., iGeppert et al.l 
l2004lPons et al.ll2009t ). 



2.3 Equation of state 

The electrical conductivity depends on the local density and 
temperature (as well as the chemical composition, which we 
take to be iron throughout for simplicity); therefore the NS 
density and temperature (radial) profiles must be known. 
We calculate the density profil e by solving the TOV stella r 
structure equations (see, e.g., IShapiro fc Teukolskvl [l983l ) . 
supplemented by an equation of state (EOS) describing the 
pressure P as a function of den sity (see, e.g.. Ilfaensel et al] 
I2OO7I : iLattimer fc Prakashll2007l. for review). We use the an- 
alytic fitTbTEaenseTT^tekhi^ (|2004l ) of the SLy EOS, 



which is a moderately-stiff EOS with a maximum mass of 
2.05 Mp) for the central d ensity pc = 2.9 x lO^'' g cm"^ 
(|Douchin fc Haensell r200lh . The NS we consider has Pc = 

I. 2 x lO^'^ gcm"^, mass M = 1.63 Mq, and radius R = 

II. 5 km; the inner crust (at p < pnuc ~ 2.8 x 10^* g cm^'') is 
at depths less than ~ 1.1 km; the outer crust (at p < pdrip ~ 
4 X 10^^ g cm~^) is at depths less than « 290 m; the heat- 
blanketing envelope (at p < pb ~ 10^" g cm~^; see below) 
is at depths less than « 100 m. The boundaries of the inner 
and outer crust and envelope are shown in Fig. [T] as well as 
the mass above a certain depth or density AM [= M — m{r), 
where m(r) is the stellar mass enclosed within radius r]. AM 
is an indication of the amount of accreted mass needed to 
sub merge the magnetic f ield down to a given density (see, 
e.g.. iGeppert et al]|l999l . for EOS-dependence of AM and 
depth). 

We only consider a single theoretical nuclear EOS and 
the standard mechanism (modified Urea process) for neu- 
tron star cooling (see below) . Varying the EOS would result 
in several effects on the Ohmic diffusion/decay timescale 
Tohm oc a{T)L^ [see eq. ([2])], in particular, a change to the 
thickness of the crust, as well as possibly inducing more 
rapid neutrino cooling for high meiss NSs. EOS effects have 
been studie d in the context of lon g-term field decay (see, e.g., 
lUrpin fc^oncnkov 1997 : Koncnkov fc Geppcrt 2001. 1 , while 
fast neutrino cool ing was needed to e xplain two of the three 
young pulsars in iMuslimov fc Pagel (|l996l ). For the CCOs, 
fast neutrino cooling is not required, and the effects of dif- 
ferent EOSs on the field evolution are beyond the scope of 
this work. 
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10 11 12 

log p (g cm-3) 

Figure 1. Relative (accreted) mass AM = M — m{r) and deptli 
i? — r as a function of density for a neutron star with tlie SLy 
EOS. Dotted lines indicate the boundary between the envelope 
and crust (at pj, Ri 10^'' g cm~^), the boundary between the outer 
and inner crust (at Pdrip ~ 4 X 10^^ g cm~^), and the boundary 



between the crust and core (at pn 



: 2.8 X 10" g cm-'^). 



2.4 Temperature evolution 

Neutrino emission is the main source of cooling during the 
first ~ 10^ yr after NS formation. The thermal history of a 
NS is primarily determined by the neutrino luminosity and 
heat capacity of the core and the composit ion (i.e., ther- 
mal conductivity) of the surface layers (se e iTsurutal 119981 : 
lYakovlev fc Pethickl [20o3 : iPage etall |2006| . for review) At 
very early times, the core cools via neutrino emission while 
the temperature of the thermally-decoupled crust remains 
nearly constant. A cooling wave travels from the core to the 
surface, bringing the NS to a relaxed, isothermal state. De- 
pending on the properti es of the crust, the relaxation time 
can t ake ~ 10 — 100 yr (|Lattimer et al.|[l993 : iGnedin et al.l 
I2OOII ). For the next ~ 10* — 10® yr, surface temperature 
changes reflect changes in the interior temperature as neu- 
trino emission continuously removes heat from the star. 

For our temperature evolution, we take the NS to cool 
by the standard (slow-cooling) modified Urea process of neu- 
trino emission, which resul ts in T ex (|Tsuruta|[l993 : 
lYakovlcv fc Pcthick 2004: Page et"alll2006l ). Specifically, we 
use (|Yakovlev et al.il201ll ) 



T{r,t) 



9.07 X lO'' K e 



2GM 



x(l + 0.12i?^)(i|^ 



,\ i/e 



(14) 



where (t>{r) is the metric function that determines the grav- 
itational redshift. For our assumed NS model (see Sec. 12. 3|) . 
the initial core temperature Tc(t = 0) = T(0, 0) = 1.1 x 
10^ K. Faster cooling processes, such as the direct Urea pro- 
cess, would cause the temperature to decrease more rapidly. 



This can lead to an increase in the electrical conductivity 
[see eq. (Illf) ]. which could slow magnetic field evolution since 
the Ohmic diffusion timescale increases. However our results 
indicate that cooling beyond modified Urea is not needed. 

Since we are concerned with NSs that are > 100 yr old, 
we assume that the NS has an isothermal core, specifically 
T{r)e'^^'^^ = constant. On the other hand, the outer layers 
(i.e., envelope) serve as a heat blanket, and there exists a 



temperature gradient from core temperatures Tc 



10" 



10 K to surface temperatures Ts ~ 10 —10 K. We combine 
the equations of hydrostatic equilibrium, 

dP_ 

dr 



(15) 



where = {1 - 2GM/c'' R)'^' ^GM/ is the surface grav- 
ity (<?s = 2 x 10^'* cm s~^ for our assumed NS model; see 
Sec. 12.31). and thermal diffusion, 



dT 
dr 



F_ 
K 



K ' 



(16) 



to obtain 



(17) 



where K is the ther mal conductivity, 
(jGudmundsson et al.lll982[ ) 

dT _ asB dP _ 3 K dP 
dp pK gfa dp 16 Qs dp ' 

where k — IGasBT"^ /3pK is the opacity due to conduc- 
tion and radiation, i.e., 1/k = + 1/kc- We consider 
the radiative opacity to be due to free-free absorption 
Kff ^ 6 X 10^^ cm^ g~^ Z^A"^pT"^/^ and electron scatter- 
ing Kcs = ZaT / ArUp. The thermal conductivity is calculated 
using CONDUCT08 (see Sec. . Note that we ignore rela- 
tivistic effects in the derivation of eq. (|17|l since uncertainties 
in the input physics exceed the effects of their inclusion on 
the results (|Gudmundsson et al.l 1 19821 ) . For simplicity, we 
assume the magnetic field does not det ermine the surface 
tempe r ature distribu t ion (s ee Sec 12. 2p . IChang fc BildstenI 
(|2004h : IChang et all (|201G| ) showed that nuclear burning 
would very rapidly remove any surface light elements for 
the high temperatures present in young NSs; therefore we 
only consider a surface composed of iron. The relation be- 
tween the interior and surfa ce temperatures is then given by 
(|Gudmundsson et al.|[l982l ) 



Ts = 8.701 X 10" K 5i 



1/4 f r[r(pb),t] 
lO^K 



11/20 



(18) 



Equation 1)17^ is solved from p = pb to lO'' g cm~^ for a 
given T\r( ph),t] from eg. (fTHl: we ta k e Pb ~ 10^° g cm~^ 
(see, e..g, iGudmundsson et al.1 1 19821 : lYakovlev fc Pethic^ 
l2004h . The evolution of the temperature profile is shown in 
Fig. [21 while the evolution of the redshifted surface temper- 
ature Ts°° Tsil - 2GM/c2i?)i/2] is shown in Fig. [3] Note 
that the temperature profile is independent of the magnetic 
field or its submergence depth or evolution. 



3 RESULTS 

Here we give examples of our solution of the induction equa- 
tion, which illustrates the evolution of a buried crustal or 
core magnetic field. Figure [2] shows the interior profiles of 
s and the normalized magnetic field br [— Bp{r)/B^,] for 
the case where the field is initially submerged at a density 
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Figure 2. Profiles of s, radial component of normalized magnetic 
field 6r, and temperature T at five time steps {t = 0, 500, 1000, 
3000, and 10* yr) for a magnetic field that is initially submerged 
at a density Paub = 10^^ g cm~^ Dashed and dotted lines are for 
the crustal and core field configurations, respectively (see text); 
the temperature profile (solid lines) is the same for both con- 
figurations. Tb (short-long-dashed line), Tm (long-dashed-dotted 
line) and Td (short-dashed-dotted line) are the magnetic (at 
B = 10^^ G), melting, and Debye temperatures, respectively. 



Figure 3. Evolution of the normalized magnetic field at the 
surface br{R) and redshifted surface temperature Tj?° for submer- 
gence densities logpi3uj3(g cm~'^) = 10.5, 11, 11.5, 12, 12.5. Dashed 
and dotted lines are for the crustal and core field configurations, 
respectively (see text); (solid line) is the same for both con- 
figurations. Triangles are the upper limits on the temperature of 
the entire neutron star surface of the CCOs Puppis A, IE 1207, 
and Kes 79, and the square is the surface temperature of the CCD 
in Cas A (see text and Table [T}. 



Psub ~ 10^^ g cm Note that p^ub — 10^^ g cm ^ implies 
an accreted mass AM ^ 3x IQ-^Mq (see Fig.[ll. For the 
crustal field configuration, one can see that the peak in the 
magnetic field at psub at t = has decreased due to dif- 
fusion/decay by about an order of magnitude in 10* yr, so 
that s and br are roughly constant as a function of density at 
p < Psub; the field also dilTuses to greater depths (p > psub) 
as the evolution continues. 

Figure [3] shows the evolution of the normalized field 
at the surface br{R) for crustal and core fields and var- 
ious submergence densities. As found in previous works 
which considered c rustal fields (|Muslimov fc Pag3 1 19951 : 
iGeppert et al.|[l999l ). the shallower the initial field is buried, 
the more quickly the surface field grows to a maximum. 
Because the peak surface field occurs at later times as the 
submergence depth increases, Ohmic field decay causes the 
value of the surface field br{R) to decrease with psub (for 
Psub ^ 10^^ g cm~"^ and t < 10* yr); thus a stronger birth 
magnetic field B, is needed to yield the same (observed) 
B. In the core field configuration, the surface field increases 
(and saturates) earlier with shallower burial. It is evident 
from Fig. [2] and [3] that the early growths of the magnetic 
field are very similar in the crustal and core field configura- 
tions. At later times, the surface field decays, and the max- 
imum field is lower in the crustal configuration than in the 
core configuration. Thus to match an observed B, a lower 
birth field B* is required in the core configuration. Finally, 
we note that the growth times shown he re are shorter than 
those seen in previous works (see, e.g., iMuslimov fc Pag3 



Il995l : iGeppert et al.|[l999h due to our use of the improved 
conductivities (see Sec. l2.2|) . and the difference increases the 
deeper the field is buried. 



4 COMPARISON TO CCOS 

4.1 Properties of observed CCOs 

We now consider the three CCOs with measured spin peri- 
ods P and either spin period derivative P or magnetic field 
(at the pole) Bp"" (see Table [l]). Note that P is obtained 
from X-ray timing, while Bp"" is determined from X-ray 
spectra, and both of these measurements contain system- 
atic uncertainties. In particular, eq. ([1} assumes a d ipolar 
magnetic field radiating into a vacuum. ISpitkovskvl (|2006[ ) 
calculated a magnetospheric-equivalent equation that yields 
Bp lower by a factor of < 1.7. On the other hand, the rela- 
tion between spectral features and the NS magnetic field is 

not definitive. 

For Puppis A. lGottheff fc HalpernI (|2009l ) find the best- 
fit model for the X-ray spectrum includes the addition of a 
Gaussian emission line at 0.79 ± 0.02 keV, though a good 
fit can also be o btained without the line. Alternatively, 
ISuleimanov et al] (|201Q ) suggest the spectral feature can be 
interpreted as an absorption line at 0.9 keV. For our pur- 
poses, we only require knowledge of the line energy, and we 
use the 0.8 keV value. For IE 1207, we interpr et the absorp- 
tion li nes seen in IE 1207 at . 7 and 1.4 keV fSanw al et al.l 
l2002l : iMereghetti et all l2002bl : iBignami et "al.. ,2003. ) as be- 
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Table 1. Central Compact Objects 





Puppis A 


IE 1207 


Kes 79 


Age (kyr) 


3.7±0.4 


~ 7" 


5.4-7.5 


P (ms) 


112 


424 


105 


P (10-1^ s s-i) 


< 35 


< 25 


0.868 ± 0.009 


Bp (10" G) 


< 4.0 


< 6.6 


0.61 




0.8-0.9 


0.7-0.8 




E (10^3 ergs s^^) 


< 10 


< 0.13 


0.30 


Tc (10^ yr) 


> 5.1 


> 27 


192 


Tg°° (lO'' K) 


< 1.7 


< 1.9 


< 3.5 


References 


1,2,3 


4,5,6 


7,8 



Notes: 

° Age estimate is uncertain by a f actor o f thre e jR oger et al.ll988(). 
Refere nces: (l) IWinkler et al.l lll98Sl'l . (2 ) ICotthelf fc Halperiil 
bOO^. (31 iGotthelf et al.l 1 120101)7^ (4) ISanwal et al. I ||2002^ , 
(5) iDe Luca et al.l | |2004|V f6TlGotthelf fc HalpernI ll2007t) . (7) 
ISun et al. I l l2004h . (8') lHalpern fc Gotthelj ||2010| '). 



ins due to the electron cyclotron resonance (jPotekIiinll2'010l : 
ISuleimanov et al.ll20ldl 'l. Spectral features due to the elec- 
tron cyclotron resonance occur at 



Ec,c = heB/m^c = 1.158 keV (B/W^^ G), 



(19) 



and we assume a gravitational redshift factor (1 — 
2GM/c^R)-^^^ in the range 1.2 - 1.35. For the two CCOs 
with spectrally-measured magnetic fields and limits on P, 
the i3p"° satisfies the upper bound set by P. 

Table [T] also shows several related parameters that are 
conventionally given for a pulsar with measured P and P: 
rotational energy loss rate 



E = 47r^/^ = 27r^7/^, 

characteristic age 
_ P_ _ 

2P " rBl' 

and braking index n 



^BpP 
Bp p 



4Bp 2 
7B3 ■ 



(20) 



(21) 



(22) 



Note that there is a selection effect of detecting sources with 
short periods since E oc P"''. The characteristic age is much 
longer than the true age if the pulsar is born with a ini- 
tial period Pi near its currently observed period P. Since 
the observed P is low for the CCOs, there is little spin pe- 
riod evolution, and 3> 10* yr. Also given in Table [T] and 
plotted in Fig. [3] are the blackbody temperatures obtained 
from XMM-Newton observations of Puppis A, IE 1207, and 
Kes 79; the tempera ture for Puppis A is an upper limit 
iGotthelf et al.ll2010l ). while the temperatures for IE 1207 
and Kes 79 a re the cooler component of the two-temperature 
spectral fits (|De Luca et al.ll2004l : lHalpern fc Gotthelj|2010l ') 
and represent upper limits on the temperature of the entire 
NS surface. Figure [3] also plots the temperature obtained 
from spectral fits to Chandra observations of the CCO in the 
Cassiopeia A (Ca s A) supernova remnant using model at - 
mosphere spectra (|Ho fc Heinkell2009l : lYakovlev et al]|201ll ): 
note that this temperature is for the entire NS surface, and 
the temperature has been seen to decrease by « 4% over the 



last 11 yr due to the cooling of the NS (|Heinke fc Hdl2010l : 
IShternin et~al]|201 j ). 



4.2 CCOs as individual sources 

We can now use our evolution calculations to constrain the 
birth magnetic field B, and submergence density psub. We 
begin by integrating eq. ([T} to obtain 



Pit) 



ft -|l/2 

J Bp[t'fdt' 



(23) 



which allows us to determine P{t). We then calculate Bp, 
E, Tc, and n. 

Figures |4]|6] show the observed values and constraints 
of the three CCOs from Table [1] The birth magnetic field 
B, is set by the normalized field fer(-R) for each submer- 
gence depth at the nominal age of the CCO (3700 yr for 
Puppis A and 7000 yr for IE 1207 and Kes 79), so that 
Bp = Bfbr{R) [see eq. Since P is directly proportional 
to Bp [see eq. ^] and E and Tc are simply related to P [see 
eqs. (|20p and (|21|l ]. the evolutionary tracks for these param- 
eters all cross at the age of the CCO. However, we see that 
Bp and n do not coincide at the age of each CCO and thus 
a measurement of either of these would allow one to distin- 
guish between the different combinations of birth magnetic 
field and submergence depth. 

For a purely crustal field and shallow submergence 
(psub ^ 3 X 10^^ g cm~^), the magnetic field decreases af- 
ter t a; 2 x 10^ yr (see also Fig. [3]). The magnetic field 
decreases more rapidly the shallower the submergence, e.g.. 
Bp ~ 10** - lO^G yr"^ for p^ub « 10* g cm'^ On the other 
hand, for pnuh ^ 10^^ g cm^^, the magnetic field increases 
more rapidly the deeper the submergence. 

For a core field and deep submergence (psub 3 x 
10^^ g cm~^), the surface magnetic field is still growing, and 
the evolutionary tracks are almost identical to the case of 
the purely crustal field. However, for shallow submergence, 
the surface field grows more quickly and saturates on the 
timescales considered here (< 10'* yr; see Fig. |3}. An exam- 
ple for the case of Puppis A is shown in Fig. [7] Note that the 
field will decay on the much longer timescale set by Ohmic 
decay in the core [see eq. ([2])]. 

Our results are summarized in Fig.[Sl where we show the 
inferred birth magnetic field B, as a function of submergence 
density psub and accreted mass AM for the three CCOs, as- 
suming the field is either confined to the crust or determined 
by the core. A measurement of both the sign and (non-zero) 
magnitude of Bp can yield a unique solution for B* and 
Psub or AM. A negative Bp can only produced in a purely 
crustal field configuration. A large positive Bp is the result 
of deep submergence in both the crustal and core fields, 
with the magnitude being the same in both cases. There is 
a small density range [psub ~ (1 — 3) x 10^^ g cm~^, where 
Bt « (1 — 3) X 10** G] in the crustal field configuration where 
Bp = 0. For the core field. Bp = for psub < 10** g cm"^, 
such that the birth field is the currently measured surface 
field [B. = (6- 9) X 10*° G]. 



© 2011 RAS, MNRAS OOO.nifTOl 



Magnetic field evolution in CCOs 7 




■D, 

W -17 



Pupisis 'a ' ^ 





3000 6000 9000 
t (yr) 



3000 6000 9000 
t (yr) 




■0, 

ra -17 



3000 6000 9000 
t (yr) 



-| T I I 1 ^ 

Puppis A 




I ^' I f-. i I 




_l 1 1 L_ 



3000 6000 9000 
t (yr) 



Figure 4. Pulsar observables: magnetic field at the pole Bp, spin 
period derivative P, rotational energy loss rate E. and character- 
istic age Tc for the CCOs Puppis A, IE 1207, and Kes 79; squares 
are values obtained from timing measurements, triangles are up- 
per/lower limits obtained from timing, and stars are from spectral 
measurements (see Table [!}• Lines are evolution models with an 
initial spin period Pi = 0.112 s and buried crustal field at density 
Psub and birth magnetic field _B,: [log Psub(g cm~'^), S» (G)] = 
(8.0,2.8 X 10^3; solid), (12.0,5.2 X lO"; long-dashed), (12.5,1.2 X 
10^2; short-dashed), and (13.0,1.5 X lO^^; dotted). Also shown 
are the calculated magnetic field derivative Bp and braking index 
n; the singularity is because of a change in sign of Bp, from a 
growing to decaying Bp. 



Figure 5. Pulsar observables: magnetic field at the pole Bp, spin 
period derivative P, rotational energy loss rate E. and character- 
istic age Tc for the CCOs Puppis A, IE 1207, and Kes 79; squares 
are values obtained from timing measurements, triangles are up- 
per /lower limits obtained from timing, and stars are from spectral 
measurements (see Table [TJ. Lines are evolution models with an 
initial spin period Pi = 0.424 s and buried crustal field at density 
Psub S'lid birth magnetic field [log Psub(g cm~'^), B» (G)] = 
(8.0,4.1 X 10^3; solid), (12.0,4.0 X lO^^; long-dashed), (12.5,4.6 X 
lO^l; short-dashed), and (13.0,5.6 X 10^^; dotted). Also shown 
are the calculated magnetic field derivative Bp and braking index 
n; the singularity is because of a change in sign of Bp, from a 
growing to decaying Bp. 



5 DISCUSSION 



4.3 CCOs as an unified source 

Taken together, the three CCOs show a non-monotonic 
trend (based on age) in their periods and magnetic fields. 
However, there is large uncertainty in the age estimates. 
Therefore let us consider the CCOs as a unified source at 
three different evolutionary epochs, e.g., the age of IE 1207 
could be « 5000 yr and the age of Kes 79 could be 
~ 6 — 7000 yr. Since field decay is required in order to 
match the observed surface fields and (assumed) ages of the 
three CCOs, core field configurations cannot provide a solu- 
tion. In Figure |9] we fit a single (crustal field) evolutionary 
track to the three measured magnetic fields. We see that for 
submergence densities psub > 3 x 10^^ g cm~'^, the tracks 
are incompatible with the fields of all three CCOs. Since 
the birth magnetic field is anti-correlated with submergence 
depth in the decaying field regime (see Fig. [S]), we obtain 
the constraint i3» > 6 x 10^^ G. Note that the maxima at 
Bp ~ 8 X 10^° G for logp^ub = 11.4 and 11.5 are produced 
from Bt = 6.2 x 10^^ G; the lower surface values are the 
result of field decay in the NS interior (see Sec. [3|. Finally, 
measuring Bp can provide a unique solution to B, and p^^h 
or AM (see Sec. [472]) . 



We have solved the induction equation to determine the 
evolution of the internal magnetic field of a neutron star, 
in order to model the field behavior (especially growth of 
the surface field by diffusion) in neutron stars with ages 
< 10"* yr. The field has a birth magnitude « B, and is 
initially submerged below the neutron star surface at a 
density Psub. We considered both a field that is confined 
purely to the crust and a field that saturates to a constant 
level in the core; the latter is valid since the Ohmic decay 
timescale in the core is ^ 10* yr. Our study builds upon 
and i mproves previous w orks (see, e.g., iMuslimov fc Pagg 
1 19951 : iGeppert et al. I ll999t ) by making use of the latest cal- 
culations of the electrical and thermal conductivities of 
matter in the de nsity-temperature reg i me relevant to neu- 
tron star crusts (Potekhin et all ll99E ; ICassisi et all l2007l : 
IChugunov fc Haenselll2007l ), and we applied the results to 
interpret observations of the recently-recognized class of cen- 
tral compact object neutron stars, which have magnetic 
fields (B « 10" G) that are lower than those detected in 
most pulsars. 

For the three CCOs with measured spin period deriva- 
tive P or polar magnetic field strength Bp, we showed that 
there is a well-defined relationship between the birth mag- 
netic field B« and submergence density paub (or accreted 
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Figure 6. Pulsar observables: magnetic field at the pole Bp, spin 
period derivative P, rotational energy loss rate E. and character- 
istic age Tc for the CCOs Puppis A, IE 1207, and Kes 79; squares 
are values obtained from timing measurements, triangles are up- 
per/lower limits obtained from timing, and stars are from spectral 
measurements (see Table [T}. Lines are evolution models with an 
initial spin period Pi = 0.105 s and buried crustal field at density 
Psub and birth magnetic field _B,: [log Psub(g cm~'^), S» (G)] = 
(8.0,3.2 X 10^3; solid), (12.0,3.2 X lO"; long-dashed), (12.5,3.7 X 
10"; short-dashed), and (13.0,4.5 X lO^^; dotted). Also shown 
are the calculated magnetic field derivative Bp and braking index 
n; the singularity is because of a change in sign of Bp, from a 
growing to decaying Bp. 

mass AM) in the case of a purely crustal field; a well-defined 
relationship also exists for the case of a core field but only at 
Psub ^ g cm~"^. Population synthesis analyses indicate 
a normal distribution for birth spin periods Pi and a lognor- 
mal distribution for birth magnetic fields: If no long-term 
field decay occurs, then the peak and width of the distri- 
butions are Pi = 300 ms and a = 150 ms and logP «(G) = 
12.95 with cr = 0.55 jFaucher-Giguere fc KaspillioOT ). while 
simulations including field decay result in Pi = 250 ms and 
ff = 100 ms and log P, = 13.25 with cr = 0.6 (|Popov et al.l 
I2OIOI ). If the CCOs are born with birth magnetic fields that 
follow these distributions (e.g., logB, > 12 at 2a), then the 
field is buried at large depths (paub ^ 2 x 10^^ g cm~^ or 
AA/ > lO~^Af0) for both the crustal and core field con- 
figurations or at shallow depths in the crustal field con- 
figuration (psub < 10" K cm~ ^ or A M < 10~^Mq). In- 
cidentally, ICotthelf fc HalpernI |20o3) argue that the slow 
spin and weak magnetic field of IE 1207 make it possible 
for accretion from a debris disc left over after the super- 
nova. However, recent optical and infrared observations of 
IE 1207 place strong limits on this disc, including an (model- 
dependent) estimate of the initial and current total disc mass 
of < IQ-e Mq and < 10"^° M©, respectively (|De Luca et all 
I2OIII ). which implies that the magnetic field could only be 
buried to a maximum density of ~ 10^^ g cm~'^. 

We showed that different combinations of birth mag- 




3000 6000 9000 3000 6000 



t (yr) t (yr) 

Figure 7. Pulsar observables: magnetic field at the pole Bp, spin 
period derivative P, rotational energy loss rate E. and character- 
istic age Tc for the CCOs Puppis A, IE 1207, and Kes 79; squares 
are values obtained from timing measurements, triangles are up- 
per /lower limits obtained from timing, and stars are from spectral 
measurements (see Table [TJ. Lines are evolution models with an 
initial spin period Pi = 0.112 s and buried core field at density 
Psub a^nd birth magnetic field P»: [log Psub(g cm~^), _B« (G)] = 
(11.0,1.2 X 10^^ solid), (12.0,2.8 X 10"; long-dashed), (12.5,1.1 X 
10^2; short-dashed), and (13.0,1.5 X lO^^; dotted). Also shown are 
the calculated magnetic field derivative Bp and braking index n. 



netic field and submergence depth can lead to the same ob- 
served values of the spin period derivative, as well as surface 
magnetic field, rate of rotational energy loss, and character- 
istic age. Predicted values of the braking index or rate of 
change of the magnetic field Bp could distinguish between 
these combinations. Although the braking index is not likely 
to be detectable in the CCOs, it may be possible (though dif- 
ficult) to measure the rate of change of the magnetic field Pp 
from the spectral lines (in the case of Kes 79, where the mag- 
netic field has only been determined from a P measurement, 
an electron cyclotron line would occur at ~ 0.5 — 0.6 keV). 
Only a crustal field configuration can produce a negative 
Pp, while a large positive Pp indicates deep submergence, 
irrespective of field configuration. If Pp is constrained to be 
very small, then the submergence density is either shallow 
for a core field or ^ 10^^ g cm'^ (AM 5 x 10"^Mq) for 
a crustal field. 

If the three CCOs are treated as a single source at differ- 
ent epochs, then the surface field is seen to decay with time, 
which rules out a core field configuration since the Ohmic de- 
cay timescale is much longer than 10* yr. For a purely crustal 
field, we showed that paub ^ 3 x lO^'^ g cm""^ and P« > 
6 X 10^^ G. A constraint at the other end can be obtained 
by considering another member of the CCOs, the neutron 
star in the Cassio peia A supernova remnant, which has an 
age of 330 ± 19 yr (jPesen et al.ll2006h . The field of the Cas A 
ceo was found to be < 10^ G from fits to its X-ray spec- 
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Figure 8. Intrinsic neutron star magnetic field as a func- 
tion of density at submergence p^ub (top panel) and accreted 
mass AM (bottom panel) for the CCOs Puppis A (solid lines), 
IE 1207 (dashed lines), and Kes 79 (dotted lines). The upper and 
lower sets of lines arc for the crustal and core field configurations, 
respectively. 



trum (|Ho fc Heinkell2009l : lHeinke fc Hdl2010l : IShternin et all 

|2011^: the no ii-dctcct ion of pu lsations from this source 
llMurrav et"al.. 2 002: Mereghetti ct al. 2002 a: _ Raiisoni 200^ ; 
IPavlov fc Lunal l2009l : lHalpern fc Gotthelj f2010l ') also indi- 
cates that the field is too low to produce pulsar-like emission. 
Its low field at a young age suggests Psub > 10^^ g cm^^ and 
-B* < 9 X 10^^ G. A better determination of the ages would 
(in)validate the CCOs as a unified source and could allow 
for a stronger constraint on their birth magnetic fields. 

We also note that, irrespective of our magnetic field evo- 
lution calculations, the P that is inferred from the spectral 
fines of Puppis A and IE 1207 resufi in _E « 4.4 x 10^^ 
and 1.7 x 10^'^ ergs s^^, respectively, while E — 3.0 x 
10'^^ ergs s^^ from the timing of Kes 79. Standard neutron 
star cooling (see Sec l2.4p yields the redshifted bolometric lu- 
minosity I/J^i, so that E/L^^i ~ 0.2 — 0.3 for Puppis A and 



Kes 79 and a much lower value of ~ 10 for IE 1207 (due 
to its slower spin period). Furthermore, E^^'^/d'^ ~ 10"^ 
(relative to Vela) for Puppis A and 10"^ for IE 1207 and 
Kes 79, where d is the source distance; these are far be- 
low the values from so urces that have been detected in the 
gamma- rays bv Fermi (fSmith et al.ll2008l : lAbdo et al.ll20ld : 
see also lZane et al.r201ll). 



Finally, iHalpern fc Gotthelj (|2010l ): iKaspi I (|2010l ) 



noted that CCOs occupy an underpopulated region in P—P. 
On the one hand, we have shown that the CCOs may just 
be representative of the low end of the distribution of birth 
magnetic fields, with B ~ 10^^ G. On the other hand, CCOs 
may have higher magnetic fields that have been submerged 
to great depths. In this case, P is increasing rapidly, and the 
CCOs are evolving to join the majority of the pulsar popu- 



Figure 9. Pulsar observables: magnetic field at the pole Bp, spin 
period derivative P, rotational energy loss rate E. and character- 
istic age Tc for the CCOs Puppis A, IE 1207, and Kes 79; squares 
are values obtained from timing measurements, triangles are up- 
per /lower limits obtained from timing, and stars are from spectral 
measurements (see Table [TJ. Lines are evolution models with an 
initial spin period Pi = 0.105 s and buried crustal field at density 
Psub a^nd birth magnetic field S«: [logpsub(g cm~'^), Bt(G)] = 
(7.6,5.1 x 10^3; solid), (11.0,9.0 x 10^^ long-dashed), (11.4,6.2 x 
10"; short-dashed), and (11.5,6.2 X 10^^ dotted). Also shown 
are the calculated magnetic field derivative Bp and braking in- 
dex n; the singularities are because of a change in sign of Bp, 
from a growing to decaying Bp. The dot-dashed line for E is the 
same as the dotted line but reduced by a factor of 4"^ due to the 
significantly longer spin period of IE 1207 [see eq. II20I I] . Note 
the different scale of the y-axes compared to those in Figs. I4I7I 
especially the linear (versus logarithmic) scale of the Bp-panel. 



lation at longer spin periods, higher P, and higher observed 
magnetic fields. 
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